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effective on-site interactions expressing Coulomb repulsion and Hund's rule. We 
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1. Introduction 

Ferromagnetism is known to originate from strongly correlated states of quantum me- 
chanical electrons with a very large total spin but small total energy. Microscopic mecha- 
nisms giving rise to a coexistence of metallic behavior with ferromagnetic ordering at suffi- 
ciently low temperatures have been proposed, but are not well understood, mathematically, 
yet. Ferromagnetism in conjunction with a positive electric conductivity (metallic behav- 
ior) is a collective phenomenon arising from a competition of spatial motion (hopping) 
of quantum mechanical charged particles with half-integer spin obeying Pauli 's exclusion 
principle, i.e., of electrons, in a periodic background with Coulomb repulsion between pairs 
of such particles. A mathematically precise description of this phenomenon is difficult be- 
cause there are two kinds of gapless excitations: electron-hole pairs very close to the Fermi 
surface, and spin waves in ferromagnetically ordered spin configurations. In a perturba- 
tive analysis of states of very low total energy, there are therefore two sources of infrared 
divergences, or of 'small energy denominators', namely electron-hole excitations with an 
energy arbitrarily close to the groundstate energy, and spin waves, or magnons, of very 
long wave length. 

In this paper we study some tight-binding models of metallic compounds with two bands 
partly filled with electrons. In a ground state, the low-lying band is at most half-filled, 
due to strong on-site Coulomb repulsion between pairs of electrons in a spin-singlet state, 
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while the higher-lying band is assumed to be weakly filled, so that a Fermi surface is 
expected. It is assumed that an electron from the low-lying band and an electron from 
the higher- lying band occupying the same site of the underlying lattice obey Hund's rule, 
i.e., their interaction energy is minimized if they form a spin-triplet state. A two-electron 
spin-triplet state is symmetric under exchange of the two spins. The Pauli principle then 
requires the microscopic orbital wave function of the two electrons to be antisymmetric 
under exchange of their positions, which makes the Coulomb repulsion between the two 
electrons small. (Concretely, an antisymmetric microscopic wave function for two electrons 
moving in the field of an ion may be constructed as a Slater determinant of, e.g., two 
different d-shell orbitals.) It may be of interest to remark that a mathematically rigorous 
derivation of Hund's rule in atomic physics from first principles has not been found, yet. 
That it is assumed to hold in our models must therefore be considered to be a reasonable, 
but heuristic ansatz. In order to eliminate small energy denominators due to spin waves 
of very long wavelength, we choose the couplings between the spins of an electron from 
the lower band and of an electron from the higher band occupying the same site to be 
anisotropic. Hund's rule cannot be invoked to justify this feature. Our results can be 
expected to hold for isotropic spin-spin couplings, in accordance with Hund's rule, too; 
but we have not been able to rigorously prove this. 

The models studied in this paper are arguably the simplest physically relevant models 
in which metallic ferromagnetism can be exhibited. 

The feature that magnetic ordering emerges from a cooperation between electron hop- 
ping and local, on-site electron-electron interactions appears to be inherent in several 
tight-binding models and has been studied intensely. Antiferromagnetic ordering, for 
example, can be seen to emerge in the half-filled (single-band) Hubbard model in second- 
order perturbation theory in t/U, where t measures the strength of hopping and U the 
strength of on-site Coulomb repulsion. This was discovered by Anderson ^ ; (for a more 
precise analysis see also A model simpler than the Hubbard model is the Falicov- 
Kimball model. There are many rigorous results about the Falicov-Kimball model starting 
with IH^J. A fairly systematic perturbative analysis can be found in 0, and references 
given there. These and many further results show how long-range correlations in ground 
states or low-temperature equilibrium states can arise from a cooperation between electron 
hopping and on-site interactions. Unfortunately, the perturbative methods in ^2 
can only be applied to the analysis of insulators, i.e., of states where electrons are essen- 
tially localized, because, in such situations, there are no small energy denominators or 
infrared divergences. The analysis of long-range correlations in metals calls for different, 
non-perturbative methods. 

One approach towards understanding simple examples of itinerant ferromagnetism is 
based on studying the Hubbard model on very special lattices that give rise to a macro- 
scopic degeneracy of the ground state energy of the Hubbard Hamiltonian (i.e., to a 'flat 
band', or to a nearly flat band); see ^l^j. A Hartree-Fock study of the Hubbard model 
also provides useful insights |21E1- More recently, there has been some interest in study- 
ing manganites described by models with several bands. Numerical studies indicate that 
there is ferromagnetic ordering at sufficiently low temperatures; see, e.g., jOnlOnMj, and 
references given there. The models studied in this paper are similar to models used to 
describe manganites. They have ground states involving macroscopically large domains of 
essentially free electrons but with aligned spins. Electrons are discouraged to leave such 
a domain, because, in neighboring domains, the direction of their spin violates Hund's 
rule, and this would result in a state of comparatively large energy. At the technical 
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level, our analysis is very much based on ^ (see also |H1E])- The methods developed 
in these references enable us to prove lower bounds on the ground state energy which, 
when combined with rather simple variational upper bounds, lead to the conclusion that 
the boundaries between domains of electrons of opposite spin have a total length growing 
much less rapidly than the total number of sites in the system, and this enables us to 
exhibit ferromagnetic ordering in the ground states. The fact that there are only two 
preferred spin orientations, ] and [, is, of course, due to the anisotropy in the on-site spin- 
spin coupling for two electrons in different bands occupying the same site. (The analysis 
of isotropic models would be considerably more difficult.) Within large domains of a fixed 
preferred spin orientation, electrons are completely delocalized, so that the ground state 
energy is not separated from the energies of excited states by a uniformly positive energy 
gap. 



Acknowledgments. We are grateful to a referee for very useful comments and sugges- 
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2. Setting, models, and summary of main results 

We consider a lattice model with electrons in two different bands, indexed by a = 1, 2. 
We let Cao-(a^) and Caa{x) denote the creation and annihilation operators for an electron of 
band a and spin a =T, i, at site x G Z*^. The state space of the system in a finite domain 
A C Z'^ is the Hilbert space 

nK = TK® Ta, (2.1) 
where Ta is the usual Fock space for electrons. 



^a = ^'(A) ® 



(2.2) 



N >0 



Here P_ is the projector onto antisymmetric functions. The energy of electrons is partly 
kinetic and partly due to interactions among themselves. The kinetic energy is repre- 
sented by standard hopping terms. Interactions are of two different kinds. First, Coulomb 
interactions are approximated by on-site operators of the Hubbard type. Second, a pair 
interaction involving the spins of electrons of different bands reminds of the Hund rule for 
the filling of atomic eigenstates. Precisely, we define the '2-band Hund-Hubbard Hamil- 
tonian' as 



H 



HH 



a=l,2 xGA 



a=l,2 cr=ta a;,j/eA 
\x-y\=l 



+ Ui2 ^ni(x) n2(x) - Si{x) ■ S2{x). 



(2.3) 



xeA 



xeA 



The parameters ta control the kinetic energy of electrons of each band and they are related 
to the effective mass of electrons. We suppose that ti > t2 and refer to electrons of band 
1 as 'light' (they move fast) and electrons of band 2 as 'heavy' (they move slowly). The 
parameters Ui,U2, U12 are positive and represent the Coulomb interaction energy between 
two particles of band 1, two particles of band 2, and one particle of each band, respectively. 
The number operators are defined by naa{x) = c}ia{x)cacr{x) and na{x) = na^{x) + nai{x). 
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Finally, J > measures the strength of the coupling between the spins; spin operators are 
given by 

S^J\x) = Y,r^£,cl{x)ca.,{x\ (2.4) 

where t'^^\ j = 1,2,3, are the Pauli matrices ^(iq), ^( ^ g'), ^(o -i)- We will consider 
a simpler model with anisotropic spin interactions (the third components of the spins 
interact more strongly) and with the property that one of the two bands is nearly flat. 
The Hamiltonian of this model is introduced in (|2.9() . below. The following discussion 
applies to general Hund-Hubbard models, isotropic and anisotropic ones. 

We note that neither the hopping terms alone nor the Hund couplings alone give rise 
to global magnetization. Without on-site interactions the ground state favored by the 
hopping terms is not magnetic; the kinetic energy is minimized by a state where half the 
electrons have spin |, and half the electrons have spin |. As for the Hund couplings, 
they are local (on-site) and do not yield the long-range correlations that are present in a 
ferromagnetic state. Ferromagnetism in this model results from a cooperation of the two 
terms. 

This model was numerically studied in [U] for ti = t2- Ferromagnetic phases were 
identified at low temperatures and for some intervals of electronic densities. The present 
model with ti < t2 actually behaves more like the two-band Hubbard model with Kondo 
spins of jni; heavy electrons here play a role similar to that of Kondo spins. Ferromagnetic 
phases were also numerically observed for this model. 

In this paper, we study a phase with spontaneous magnetization. The total spin oper- 
ator in a domain A C Z'^ is denoted Ma and is given by 

Ma = Y,[Si{x) + S2{x)]. (2.5) 



The norm of Ma is 



Ma|| = (E(<)^)'^'> (2-6) 



(i) 

with Af^ the i-th component of Ma- We expect that the system displays extensive 
magnetization. That is, in a domain A the expectation of H-ZWaII in the ground state 
should be proportional to |A|. We are not able to prove this, but we can prove that the 
system is magnetized at least on a 'mesoscopic' scale. To be precise, we consider the 
following definition of the magnetization per site: Given a state T E TCa and a subdomain 
A' C A, we define 

"^^.^' = m Yl r^(T^'II^A'+x||T); 

x:A'+xCA ^2 7) 

' ' x:A'+xCA ' ' 

Note that we consider the expectation value of the norm of Ma' in the state T, averaged 
over all translates of A' such that A' + x remains in A. The number of terms in the sum is 

(3) 

|A|, up to a negligible boundary correction. We clearly have that tua^a' ^ ""^a a" ^^^^ 
prove that, for the anisotropic model introduced in (|2.9|) . below, the third component of 
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the ground state magnetization per site, rn^l^r, satisfies 

(3) 

lim lim m\ . , > 0. 

A' /'Z'i A yzd ' 

The hmits are over boxes of increasing size, and electron densities are kept constant. 

Perturbation methods shed some light onto the structure of the phases of this model. 
The situation is similar to the one in the Hubbard model at half-filling and with strong on- 
site repulsion, which behaves like the antiferromagnetic Heisenberg model. Perturbative 
methods along the lines of [HJ [3] can be used for many rational densities. However, such 
methods can be applied only when electrons are localized, and this is not the case in a 
conducting metal. Here we focus our attention on regimes where there does not exist an 
energy gap separating excitations from the ground state energy, and where some electrons 
have delocalized wave functions. 

The Hamiltonian (|2.3|1 is symmetric with respect to spin rotations, and this makes the 

(3) 

study difficult. We simplify the model by replacing the term Si{x) ■ S2{x) by S{ (x) ■ 

(3) 

S2 (x). Furthermore, we let t2 — > 0. We can fix the energy scale by choosing ti = 1. We 
then obtain the following simplified Hamiltonian ('2-band Ising-Hubbard model') 

nai{x)nai{x) + Ui2'^ni{x)n2{x) 

o" x,y£A a=l,2 xGA a;GA 

l^-y|=i (2.8) 
~ 4 E[^iT(a^) - nilix)] ■ [n2t(x) - n2iix)]. 

zGA 

Let pi, p2 denote the densities of light and heavy particles, respectively. We prove 
the following statement in Section^ Our proof works in dimension larger or equal to 2 
(Equation (|3.19|) holds for d ^ 2 only). We do not have results for the one-dimensional 
models. 

Theorem 2.1. Let d ^ 2. For arbitrary < pi < p2 ^ 1, there exists Jq < 00 such that, 
for J — U12 > Jo, all ground states T of Hj^ satisfy 

lim lim m^^\, = ELltJ^_ 

This theorem suggests that the ground state displays 'saturated ferromagnetism', as 
it has maximum total spin. Notice that for large enough J it holds independently of 
Ui,U2^0. 

The proof of Theorem 12.11 reduces to the study of the ground state energy for fixed 
configurations of heavy electrons, since the latter do not have kinetic energy. To a con- 
figuration of heavy electrons we can assign domains of | and J, spins. A light electron of 
spin I, say, is in a state that is essentially localized on the domain where heavy electrons 
also have spin |. Hund interactions suppress other configurations. In the limit J ^ 00, 
the ground state energy is purely kinetic and it is minimal in a configuration of heavy 
electrons with large domains of identical spins. This allows us to show that the size of the 
boundary of these domains is less than |A|^~d , meaning that spins are locally aligned. See 
Section 0] for details. 

The constant Jq in Theorem 12 . II depends on pi, and there are good reasons for it. There 
is no ferromagnetism for pi = or pi = p2 = 1. The former case results in independent 
spins at each site. The latter case can be treated with perturbation methods. Non-empty 
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sites are typically occupied by two particles of spin |, or by two particles of spin |. An 
effective interaction of strength --^-r stabilizes antiferromagnetic chessboard phases in 

(7l+2 

the ground state and at low temperatures. This interaction can be obtained using the 
method described in jHI- The case pi = p2 < 1 is more subtle. 

We can improve the result of Theorem 12.11 and consider a model that interpolates 
between ()2.3p and (|2.8|) . We refer to the following Hamiltonian as 'the asymmetric 2-band 
Hund-Hubbard model': 

^A™ =-T.taYl E cL(^)Caa(y) + Yl UaY^^ari^Mx) 
a=l,2 (T=t,J. x,yeA a=l,2 xGA 

\x-y\=l 

+ Uu "1 (^) ^2 ix)-jY (^) (2.9) 

zgA xeA 

xGA 

Let us again fix the energy scale by setting ti = 1. 

Theorem 2.2. Let d ^ 2. For arbitrary < pi < p2 ^ 1, there are constants Jq < oo and 
c > (both depend on the densities) such that if ^ — U12 > Jo o,nd t2, J"*" < c, all ground 
states T of Hf^^ satisfy 

lim hm mAA, = ^l±^. 
K' /TA A ' 2 

This theorem again holds uniformly in Ui,U2 ^ 0. 

Our paper is organized as follows. We discuss the properties of the ground state of a 
simple model in Section |31 The results for the simple model are then used in Section |^ 
where Theorem 12.11 is proved. Finally, it is shown in Section [3 that the claims for the 
Ising-Hubbard Hamiltonian (|2.8j) can be extended to certain perturbations, that include 
the asymmetric Hund-Hubbard model (|2.9)) . This proves Theorem 12.21 

3. Interacting electrons in a magnetic potential 

We introduce in this section a Hubbard model of electrons in an external potential that 
involves the third components of the spins. We do not insist on the physical relevance of 
this model. The sole motivation for this section stems from applications to Hund-Hubbard 
systems. We will use Propositions 13 . 1^3 . 31 in Sections |1] and [5] in order to prove Theorems 
12.11 and l2. 21 — these theorems being physically motivated. 

The results below extend the bounds for the ground state energy of spinless electrons 
in binary potentials proposed in "T. We work in the Fock space J-\ of spin ^ electrons 
in A. Let F be a 'magnetic potential', that is, y is a collection of nonnegative numbers 
Vx , Vx indexed by sites x G Z . The Hamiltonian is 

HAiV) = - Y E ct(xK(y) + C/j;nt(xK(x) + ^ K^n.(x). (3.1) 

Here, ct{x) and Ca-(x) are creation and annihilation operators of fermions of spin a at x, 
and no-(x) = ct(x)co-(x). We suppose that a gap Vq separates the minimum value from 
other values of the potential. Introducing 

A„ = {xeA:V^ = 0}, A = Ai^uA^ (3.2) 



x,y€A cr=t,i xGA a:GAcr=T,i 

\x-y\=l 
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(the sites where the potential is zero for some spin), we define 

Vo = min inf . (3.3) 

We assume that Vq is strictly positive. 

In order to understand the bounds on the ground state energy given below, it is useful 
to consider the situation where Vq — > oo. Assuming that j4-|- n = 0, the domain A is 
partitioned into and A \ A. Electrons of spin a are described by wave functions 

with support in the energy being infinite otherwise. Electrons do not interact and 
their ground state energy is purely kinetic. It mainly consists of a bulk term that depends 
on the electronic density inside A^- and that is proportional to the volume \Aa\. The effect 
of the boundary of A^j is to increase the ground state energy by a term proportional to 
the size of the boundary. The ground state energy of non-interacting spinless electrons in 
arbitrary finite domains was studied in [7j; upper and lower bounds were established that 
confirm the discussion above. As Vq decreases from infinity to a finite value, electrons 
delocalize somewhat, but the situation does not change in any essential way. 

Estimates for the ground state energy involve the energy density of free spinless electrons 
in the limit of infinite volume. As is well-known, the energy per site e(p) for a density 
< p < 1 of electrons is given by 

^(P) = -prXd £kdk, Ek = -2^cos ki, (3.4) 

l^^i j£fe<£F(p) -^^ 

where Epip) is the Fermi energy, defined by the equation 

P=7^f d/c. (3.5) 

Notice that e{p) < for < p < 1. We need to define the boundary B{A) of a set ^4 C Z*^; 
it is convenient to define it as the number of bonds that connect A with its complement, 

B{A) = #{{x,y):x£A,y(^ A, \x - y\ = 1}. (3.6) 

We first give a bound for fixed densities of electrons of each spin. In the absence of 
interactions {U = 0) the following proposition merely rephrases similar results in [2]. We 
define E\{V;N-^,N^) as the ground state energy of H/^{V) when the number of spin | 
(spin I) electrons is N-^ {N^ respectively). We introduce a notation for electronic densities 
inside A^ and A^; for a =t,i, we let p^j = We have the following bounds for the 

ground state energy. 

Proposition 3.1. Let V be a magnetic potential and N-^ , Ni be numbers with the prop- 
erties that: 

• A^nAi = $. 

• Vo> 2d{Vd + l). 

• iV| ^ \A^\, ^ \A^\. 

Then there exists a{p) > (independent o/Vq), for < p < j^, such that 



[eip.)\A,\ + {aip„) - ^iVo))BiA,] 
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The inequalities in this proposition hold uniformly in U. The proof of Proposition 
is based on results in ref. [Zj, where the sum, Sa,n, of the N lowest eigenvalues of the 
discrete Laplacian txy = —S\x-y\,i,x,y G A is estimated, with A a finite set of lattice 
points of arbitrary shape. Two of the results in jjj are relevant for our analysis: 

• We have upper and lower bounds, 

e(p)|A| - '-^B{A) ^ Sa,n ^ e{p)\A\ + a{p)B{A), (3.7) 

where p = j^, and a{p) is strictly positive for any < p < 1. Recall that e{p) 

is negative, so that all boundary terms in the above equation are positive. (The 
notation in is slightly different, the Hamiltonian being shifted by 2d and the 
boundary is defined differently.) 

• If S^j^ denotes the sum of the lowest eigenvalues of the operator —S\x-y\^i + 
UX\c(^x), where Xac is the characteristic function of the complement, A'^, of the 
set A, and if U is positive, we have 

Sa,n ^ S% ^ Sa,n - 7{U)B{A), (3.8) 

for some ^(U) — > as C/ ^ oo. 

The upper bound for E\{V; N^^, Ni) does not depend on Vq. Increasing the values of the 
potential actually incresases the energy, so it is enough to prove the statement in the limit 
Vb ^ oo. Electrons with different spins are independent, and the upper bound follows 
from the one in (|3.7|) . 

Let us turn to the lower bound. The operator that represents interactions between 
electrons is positive; we get a lower bound for the ground state energy by taking f/ — > 0. For 
Vq = oo we are in the situation of [7]. For finite Vq we use ()3.8|) with minor modifications. 
Namely, starting with Equations (4.3)-(4.5) of [7j but introducing our measure B{A) of 
the boundary, the upper bound in Eq. (4.8) can be replaced by ^^S(A). It is useful to 
modify the bound for the number of sites at distance n from the domain A (recall that we 
are using the distance here). It is not hard to check that 

#{:r:dist(x,A) = l}^S(A); 

#{x : dist (x, A) =n} s^d- #{x : dist (x. A) = n - 1} ff n ^ 2. ^ ' ' 

We therefore have that #{x : dist (x. A) = n} ^ d^'^^B(A). This allows to bound N — Tip 
in Eq. (4.11) by (^Y^_2dy2_4^^i B{A), leading to the present definition of 7(Vb). The bound 
given here is better for large Vq than the one in [7j. 

The considerations above show that Proposition 13.11 is a mild extension of |7j . The 
following proposition needs, however, a more detailed proof. 

Proposition 3.2. Under the same hypotheses as in Promsition \S.lV we have that, for all 
normalized ground states T of H\(V), 

Y^(T,n^ix)T)^^[BiA^)+B{AO]. 

(7=1,1 X^Aa 

Proof. Let M„ be the number of electrons of spin a that are outside of Ao-, and Pm^m^ be 
the projector onto the subspace spanned by states with exactly particles outside A„, 
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a =1, [. A state T can be decomposed as 

Tf" = ^ ^ CAf^M|TrM|Afp (3.10) 
Af|=0 Af|=0 

with CM|M| = ||-PM|Af^T|| ^ 0, I]c|f^A^| = 1' s-^d Tm|Mj^ = c'^j^m^Pm^m^'^ is normalized. 
The goal is to estimate 

^ ^(T,n.(x)T)= ^ 4,^^,^(M^ + Mi). (3.11) 

<T=T,ia:^^<, M|,Af| 

The strategy is to obtain a lower bound for Ea{V; N^^, Ni) that involves the expression 
above. Comparison with the upper bound of Proposition 13. II wih prove the claim. 

The ground state energy is increasing in U so that we can again set U = when 
discussing a lower bound. The Hamiltonian H\{V) can be spht into 

^a(^) = E + - E E [cUxKiy) + cUyK{x)]. (3.12) 

\x~y\=l 

Hamiltonians H^iV) consist in kinetic terms for particles of spin a in the corresponding 
domains, and of the potentials given by . They leave the subspace with fixed M| and 
M| invariant. The norm of the last operator is smaller than 2B(A^) + 2B{A^^). Therefore 

{T,H^{V)T)^ ^ cl^M,Y.('^M,M,,[HZiV)+Hl^^^{V)]TM,M,) 

-2B{A^)-2B{A^). (3.13) 

Inserting the lower bound for the sum of the lowest eigenvalues of the discrete Laplacian in 
a finite domain, neglecting the positive boundary correction term, we get the lower bound 

{Tm,m,,H^AV)Tm,m,) ^ e{^^^^)\A,\ ^ e{p„)\A,\- e^{p,)M,. (3.14) 

The second inequality holds because e{p + r]) ^ e{p) +r]eF{p) (indeed, e{p) is convex and 
its derivative is sp{p)). 

The Hamiltonian H'^^^ (V) is the second-quantized version of a one-body Hamiltonian, 
whose eigenvalues are bigger than Vq — 2d. Since £f{p) ^ 2(i, we have 

{T,HAiV)T)^ eip^)\A^\ + {Vo-M) cl^,A^,(Mt + M^) - 2i?(^T) " 2i?(^i). 

<T=Ta A/|,A/| 

(3.15) 

The right side must be less than the upper bound for E/^iV; N-^,N^) stated in Proposition 
13.11 Using — 4^ ^ 1, we get Proposition 13.21 □ 

We turn to the situation where the total number of electrons is specified, but not their 
spins. Let E\{V; N) be the ground state energy of H\{V) with electrons. In the proof 
of the following proposition we have to assume that the dimension of the system is at least 
2. 

Proposition 3.3. We suppose d ^ 2. Let V be a magnetic potential and N be a number, 
that satisfy 

• A^nAi = $. 

• Vo> 2d{Vd + l). 



10 



JURG FROHLICH, DANIEL UELTSCHI 



• iV ^ |^|. 

Let p = < 1; then there exists a{p) > such that 

e{p)\A\ - + BiAO] > Ej,{V- N) ^ e{p)\A\ + (a(p) - ^(Vo)) + BiA^]. 

Proof. The upper bound follows from the upper bound of Proposition Kill that holds for 
all U. We can set U = for the lower bound. Because electrons of different spins 
do not interact, the ground state energy is given by a sum of lowest eigenvalues of the 
corresponding one-body Hamiltonians for particles of given spin. Let N'^ be the number 
of spin I electrons in the ground state. Taking into account multiplicities, there are 
1^1^ I available eigenvalues in {—2d, 2d) for spin | electrons, and |^|| eigenvalues for spin 
I electrons. Other eigenvalues are larger than Vq — 2d. Since N ^ \A\, we must have 
^ A^l ^ l^ll and ^ iV - iV| ^ 

Let us introduce p = jj-^, p' = j^y , and rj = -^j^. We have = ^ and pi = Using 
the lower bound of Proposition 13.11 we obtain 

Ea{V; N) ^ {r?e(^) + (1 - v)e{^)}\A\ + [a(^) - 7{Vo)]B{A^) 

+ K^)-7(^o)]5(^i). (3.16) 

This bound does not hold for all p' , but it holds when p' corresponds to a ground state. We 
get a lower bound by minimizing over p' . A difficulty arises, namely that the coefficient 
of B{A-^) or of B(A^) could be negative. The term in braces reaches its minimum for 

— = = p. Let e be such that the minimizer for the whole right side be — = p — - 

(and = P + T~^)- The fractions ^ and are small, because a and 7 are small. 
Hence e is small, too. Let /(e) denote the term in braces, 

/(e) = iie{p - f ) + (1 - ^)e{p + j^). (3.17) 

The second derivative is 

f"{e) = '-s'Ap - f ) + jki4{P + T^)- (3.18) 
One easily verifies that £p{p) ^ c, with c strictly positive when d ^ 2. This implies that 

/(e) > <P) + 2^^^- (3-19) 

Then 

Ea{V- N) ^ e{p)\A\ + ^1^1 + [a{p - ^) - im]B{A^) 

+ [a{p + jh^)-^{Vo)]B{A^). (3.20) 

The right side should be e(/?)|A| + a[B{A^) + B{Ai)\. We must show that the brackets 
are strictly positive, depending on p and Vq, but uniformly in rj. Four situations need to 
be carefully investigated: (1) if p is small and | > is of the order of (2) if p is small 
and < is of the order of p; (3) if p is close to 1 and ^ < is of the order of 1 — p; 
and (4) if p is close to 1 and > is of the order of p. These four cases are similar, 
so it is enough to consider case (1). The factor in front of -B(Aj) is bounded away from 
because p + is bounded away from and 1 (uniformly in ry), so that a{p + jz") > 
(uniformly in rf). We can assume that - > f (the bound is uniform in rj otherwise), and 
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we consider the factor in front of B(A-^). We take advantage of the second term in |^|, 
observing that 

> T^(f)M^I- (3-21) 

Now r]\A\ = ^ ^B{A'^), and we see that the factor in front of is uniformly 

bounded away from zero as — > 0. □ 

4. The Ising approximation 

In this section we prove Theorem l2.1l Heavy electrons are static and they can be treated 
classically. Their state is represented by a classical spin configuration sa G {0, |,|,2}'^, 
and the model (|2.8|) corresponds to a Hamiltonian, H/<^{s/<^), acting on J^\. The expression 
for Ha{sa) is given by 1)2. 8|) . with the understanding that the operators c\^{x), cio-(x) act 
on J^A (instead of ./^a)) and the operators n2aix) are replaced by numbers as follows: 



n2-f(x) 



7121 (x) 



1 if Sa; =t 

ifSa;=i, 

'o ifs^=T 

1 if =i . 



Thus H/<^{s/<^) is a Hubbard Hamiltonian with an external potential (or 'field') given by 
SA- It is convenient to add a constant j — Uu to the energy so that the potential is 
non-negative. For given sa we define the potential by 



X 



f if =T 

i if =i 

i-Ui2 ifs, = 

' + Ui2 if s,^ = 2. 



(4.1) 



V 4 

Next, we define in the same way, by flipping the spins. With A^i = />i|A|, the Hamilton- 
ian for the Ising-Hubbard model can be expressed using the Hamiltonian H\(y) defined 
in 1)3. namely 

^a(sa) + (i - Ui2)Ni = Ha{V) + U2J2 '^-.,2. (4.2) 

The strategy of our proof of Theorem 12. H is as follows: 

• A state where all electrons have spin | gives us an upper bound for the ground 
state energy (Equation 1)4. 3|l '). 

• We derive a lower bound for the ground state energy that involves 'classical ex- 
citations' of Sa — regions where heavy particles do not have parallel spins. See 
Proposition 14.11 

• By combining the upper and lower bounds for the ground state energy, we find 
that any ground state configuration necessarily has only few excitations (see Eq. 
()4.5|) ). This suffices to prove Theorem 12.11 

Let E{s\;pi) denote the ground state energy of 1)4. 2() . and let be the number of 
heavy electrons in sa. A candidate for the ground state is a purely ferromagnetic state, 
where all particles have spin |. Heavy electrons occupy a domain A-^ with |74-[-| = A^2) and 
they are described by the configuration Sx =t for all x € A-^. Light electrons also have 
spin t; electrons of identical spins do not interact, so that the ground state is given by the 
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Pi|A| lowest eigenstates of the hopping matrix in A. By the upper bound in Proposition 
13.11 we have that 

min^(sA;pi) ^ e{^)N2+AdNl~~^. (4.3) 



The ratio ^ = ^ represents the effective density of hght electrons when they all reside in 
A-^. The second term on the right side is an upper bound for the boundary contribution 
to the energy of an optimal domain with A''2 sites. 

Next, we turn to a lower bound. It is useful to introduce 

i{p)=pE^{p)-e{p). (4.4) 

Notice that ^ ^(p) ^ 2d, and ^'(p) = p^pip) > so that ^(p) is increasing. Recall the 
definition ()3.2() for sets A-^ and A^ that are determined by the potential ()4.1|) . Notice that 
A is the set of sites occupied by exactly one heavy particle, and that A-^ Ci A^ = 0. 

Proposition 4.1. For < pi < p2 ^ 1; there exist Jq < oo and a > (both independent 
of A) such that if j — U12 > Jo, we have that 

E{sa; pi) > e{^)N2 + U^) + iC/2] {N^ - \A\) + a[B{A^) + B{A{)\ , 

for arbitrary s\. 

Remark: A similar bound can be proven when 0<pi<2 — p2^1- 
One main consequence of Proposition 14.11 can be obtained by combining it with the 
upper bound (|4.3|) . We get 

f 1^1 ^ iV2(l - const • N~^^'^) 

{ i-i (4-5) 

[b{A^) + B {A i) const ■ N2 " 

for constants that are uniform in the size of the system. These inequalities imply Theorem 
12. H as is shown below. 

Proof of Proposition We observe that '^^.Ss^,2 = ^(-^2 — |^|) yielding the term in- 
volving U2- To alleviate our notation we suppose now that U2 = 0. Let ^ > be a small 
number; we first consider configurations such that ^ 1 — ^. Proposition 13.31 gives 

EisA-, Pi) ^ e(^)| A| + (a(^) - 7(i - C/12)) \B{A^) + B{A^)\. (4.6) 

The function J^e(^) is convex in z^, and its derivative with respect to v is equal to — 
Therefore 

e(^)|^| = e(£^^)^7V2 ^ e(^)iV2 + i{^){N2 - \A\). (4.7) 

Since < 1 — C the function a(j^) is uniformly bounded away from zero, and we 
obtain a strictly positive 5, provided \ — U\2 is large enough (see Proposition 13.11) . 

We now consider configurations such that 1 — < ^ ^ 1. Equation (|4.H|) is still valid 
but a(-^) may be very small and we ignore it; it is positive. Convexity of e(p) yields 

e(f^)^e(l) + (^-l)£F(l)^ -2dC. (4.8) 

We have used that e(l) = and eF(l) = 2d. Because B{A{) + B{A{) ^ 2d\A\, we obtain 
from 

£;(sa;pi)^ - 2d[C + 7(1-^12) + 5] 1^1 +5[5(^t)+^(A)]- (4-9) 
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In order to complete the proof of Proposition 14.11 we need to check that 

-2d[C + - f/12) + 5] ^ e(£i)iV2 + e(^)(iV2 - |^|). (4.10) 

We have that N2 - |A| = A^2[l - = A^2[l - ^ + 0(C)]- As \A\ ^ N2 and because the 

term in brackets can be arbitrary small (depending on pi, /32), it is enough to check that 

0>e(^) + e(^)(l-^). (4.11) 

Using the definition (|4.4j) of ^, the condition can be reduced to ^(p) — e-p{p) > for 
< p < 1. This is easy to verify, as this function is strictly decreasing and ^(1)— eF(l) = 0. 

Finally, the case where A'^i > |^| is easy because we can use Proposition 14.11 for the 
lowest \A\ eigenvalues, and remaining eigenvalues are larger than j — U12 — 2d ^ 0. □ 

Proof of Theorem \2.1\ We have established inequalities H4.5|) that show that ground state 
configurations of heavy electrons consist of large domains with one particle of spin | at each 
site, large domains with one particle of spin |, or domains void of particles. Boundaries of 
these domains are 'sparse'. Recall that definition 1)2. 7|) of the magnetization mA,A' involves 
an average over translates of A'. It is enough to restrict to boxes that are fully in or 
in A^. Indeed, few boxes are intersecting their boundaries, and there are virtually no 
electrons outside of A. 

It is clear that rn'^^, ^ £l±£2 £qj. ^^jj states with densities pi and p2 of light and heavy 

(3) 

electrons, so that it suffices to establish the converse inequality. The definition of nvi^^, 
involves a sum over translates of A' that are inside A. All terms are positive, so we get a 
lower bound by restricting the sum to translates that are contained in either Ai^ or Ay: 

\{T,Mf\j)\^ E (4.12) 

xiA'+xCA 2;:A'+xcA| x-A'+xcA^ 



(3) 

Let us recall the definition of M^, : 



mJ^) = ^ M(3)(2;), (4.13) 



r(3) _ 

xeA' 

with 

M^'H^) = \ E (^«t(^) - nai{x)). (4.14) 

a=l,2 

Since ||M(3)(x)|| 1, we have from (jlH^ 

Y (T,|Mi^^JT)^ ^(T,M(3)(x)T)- 5](T,M(3)(:,)T)-|A'|[S(ylT) + i?(^i)]. 

x:A'+x<zK a;GA| a;GA| 

(4.15) 

Let Niy = I]a;gA(^''^iT(^)'^); then 

^(T,M(3)(x)T) = i|At| + iiVit- E(T^'^it(^)'^)- (4-16) 

xeA| x^A-^ 

The latter term is less than - — — —B{A'^) by Proposition 13.21 The same argument 

applies to spin J, electrons. Using inequalities (|4.5|) . we see that mA,A' is larger than 
\p2 + \pii up to a term of order lAI"^/*^ (it depends on A'). This term vanishes in the 
limit A ^ Z'^. □ 
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5. The asymmetric Hund-Hubbard model 

We now turn to the proof of Theorem 12 .21 The asymmetric Hund-Hubbard model (|2.9jl 
can be expressed as a perturbation of the Ising-Hubbard model (|2.8() . Namely, with ti = 1, 

Hf^ = Hj^-t2Yl ^ 4(x)c2.(y)-J^j;[5j')(x)s(^)(x)+5P(x)sf (x)]. (5.1) 

cr=t,i x,y£A xSA 
\x-y\ = l 

In the previous section we showed that any ground state configuration of the Ising-Hubbard 
model satisfies inequalities ()4.5() . Heavy electrons are now quantum particles and a classi- 
cal configuration cannot be an eigenstate. We can extend H4.5() by expanding the ground 
state in the basis of configurations of heavy particles, and show that H4.5|) holds in aver- 
age. Namely, we denote by $(sa) G the normalized state of heavy electrons in the 
configuration s\. Clearly, ($(sa)) is a basis of J^a. Any state T G ® has a unique 
decomposition as 

r = Y, c(sa)^'(sa) ^{sa), (5.2) 

where c(sa) ^ satisfies Y2sa '^^(^a) = 1, and \I'(sa) is some normalized state that repre- 
sents the light particles. Notice the asymmetry in notation: ^'(sa) is indexed by sa, but 
the configuration of spins of light particles may be very different from the configuration 
SA, in general. In particular, ^'(sa) describes a state with A'^i particles, while <&(sa) has 
N2 particles. Let X{sa) denote the number of 'excitations' of sa; namely, 

^(^a) = (iV2 - 1^1) + B{A^) + B{AO. (5.3) 

The extension of (|4.5() is as follows. 

Proposition 5.1. Let d ^ 2, and < pi < p2 ^ 1- There are constants Jq < 00 and 
7 > such that if ^ — U12 > Jq and t2,J'^ < 7, and if c{s\) are the coefficients defined 
in H5.2p for a ground state of , we have 

^c2(sa)X(sa) ^ const • A^2~^ 

SA 

for a constant that is independent of A. 

Proof. We again define V by (|4.H) and H\{V) by (|4.2|) . U2 plays no role here, just as in 
Section \^ so we set it to from now on. By the variational principle we find an upper 
bound for the ground state energy by considering a state where all electrons have spin |. 
Heavy electrons are packed together and light ones are in appropriate delocalized wave 
functions with support on A^. Neither the t2 term nor the term contributes to the 
energy of this state, and (|4.3|) therefore continues to be an upper bound for the ground 
state energy. 

The goal is now to find a lower bound with the same bulk term as in the equation above, 
plus a correction that involves the average of X{sa)- For T expanded as in 1)5. 2|) . we have 
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SA 

E c(.A)c(4)(^(^A),^'(4))(^(^A),4wc2.(y)^(4; 

SA,s'jyO-=1,l x,y£A 
\x-y\=l 

SA,SAa;eA 



(5.4) 

The first term involves the same H\{V) that appears in (|4.2)) : this gives us the right 
bulk contribution. The other two terms are actually irrelevant and it is enough to find 
estimates. 

We observe that the second term on the right side of (|5.4|) is less than 

*2 ^'c(.A)c(4)^t2(E'^'M'^'(E'^'(^A))'^'=*2 E'^'(^a)- (5.5) 

Primed sums are over configurations that are identical except for a heavy electron moved 
to a neighboring site. Given sa, there are less than 2d{N2 — \A\) + B{A^) + B{A^) such 
configurations (recall that ^{N2 — \A\) is the number of sites that are occupied by two 
heavy electrons). It follows that ()5.5() is smaller than 2(it2 X^^^ c^('Sa)-'^('Sa) (and it is 
larger than the negative of this expression) . 

The third term of the right side of (|5.4|) can be treated in the same spirit. It is necessary 
to cast the perpendicular Hund interactions in a form that shows that their contribution 
is no more than the boundary between domains of identical spins. We therefore introduce 
standard operators si~^\x), si \x), by 

Si+\x) = Si'\x) + iSi'\x), 

Sj-\x) = Si'Hx)-iSj^\x). 

Perpendicular spin interactions become 

S['\x)si'\x) + S['\x)si'\x) = 1 [s[-^\x)si-\x) + s[-\x)si-'\x)] . (5.7) 

Let X G A, and be the configuration obtained from sa by flipping the spin at x. The 
third term of (|5.4j) is equal to 

-^•^^ E E <«a)<4) (*(«a), s[*^ {x)^{sl)) , (5.8) 



SA xGA 

with ^ ='+' if =1, and # ='— ' if =T- Since s[~^\x) = c^(x)ci|(x), the Schwarz 
inequality yields the bound 

vI/(sA),5;+)(x)vI/(sX))| ^ (vI'(sA),ni^(x)^(sA))'^'(^(4),^u(^)^(4))'^'- (5-9) 

A similar inequality holds when S'|^'*(x) is replaced with s[ \x); one should simply in- 
terchange ni|(x) and ni^{x) on the right side. Then echoing (|5.5p . the absolute value of 
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(|5.8|) is found to be smaller than 



1/2 



1/2 



= 4 E E ^'(«a) [^{sK),n,,-s. (x)^(sa)) . (5.10) 

This expression is reminiscent of the expression in Proposition 13.21 However, ^(sa) is not 
a ground state for light electrons in the magnetic potential given by sa and therefore the 
proposition does not directly apply. 

We can nevertheless recycle the ideas underlying the proof of Proposition 13.21 
Let N = (A''|, A^j^, M|, Mj^) be four positive integers such that + = A^i, and 
M(j ^ N(j. No- is the number of light electrons of spin cr, and is the number of light 
electrons of spin a that are not localized on the favorable sites A„. We can expand ^'(sa) 
according to AA, in a fashion that is reminiscent of (|3.1flj) . 

c(sa)^'(sa) = E c.^(«a)^'^(sa) (5.11) 

where coefficients are positive and states are normalized. With this notation we observe 
that (|5.1fl|l is bounded above by the following expression similar to (|3.11j) 



EE^^("A)(Mt + Mi). (5.12) 



2 



The Hamiltonian H/^(y) can be split as in (|3.12jl and we obtain the lower bound (|3.13|) 
with c(sa)^('Sa) in lieu of T, and cj^{s\)'i!j\f{s\) in lieu of cm^m^'^ M^My We then get 
(|Xn|) and (|XT^ . Exphcitly, the lower bound for ifOjl is 

<^)^2 + EE^^("A){(i-C/i2-4d-4)(M| + Mi)-2S(^^)-2S(^^)-2dt2X(5A)}. 

The bulk term e{^)N2 above comes from ()3.15() . minimizing over N-^ and N^. This 
expression is less than the upper bound (|4.3|) : this implies that 

Y,Y.-"^{s^m+Mo ^ J 'l^'tf%^ E^^(^A)-^(^A)- (5-13) 

SA M I-U12-M-— SA 

(We used that N2 ^ X(sa).) Notice that the factor on the right side is small. This 
estimate is necessary to bound (|5.12|) . 

Using Proposition 14.11 for the first term in 1)5. 4(1 . we then conclude that 

Y,<?{sk)(^^{sk),Hk{V)^{sk)) ^e(£^)7V2 + 5Ec'(^A)^(*A)- (^'l^) 

SA SA 

(We assume here that 5 is smaller than ^(^) + if/g.) A gain invoking the upper bound 
(|4.3|) . we obtain 

e(^)N, + MnI--^ ^ e{fjN2+ [s - 2dt2 - ^ /l^^J^.fV ] E ^'(^a)X(.a). 

(5.15) 
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The quantity in brackets is strictly positive when j — U12 is large enough, and this proves 
Proposition 15.11 □ 

Proof of Theorem \2.<A The proof is similar to that of Theorem 12.11 except that we use 
Proposition 15.11 instead of the inequalities ()4.5() . All equations until 1)4. 14() hold without 
change. Eqs 1)4. 15() and 1)4. 16() need to be modified because the Ao-'s are not fixed here. 
These equations hold when averaged over s\ with weights c^(sa)- We obtain 

(TjMi^).JT)^ j;j;c^(sA)j;[i|A.| + liV<,-iM.-|A'|B(^.)" 

a;,A'+xCA Af a 

^ \{Ni + N2) - 2|A'| c\sa)X{sa). (5.16) 

SA 

We have used 1)5. 13() in order to estimate the contribution of Af| and M^. After division 
by |A|, the last term vanishes in the thermodynamic limit by Proposition 15.11 □ 
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